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A primordial cosmological magnetic field induces Faraday rotation of the cosmic microwave back- 
ground polarization. This rotation produces a curl-type polarization component even when the 
unrotated polarization possesses only gradient-type polarization, as expected from scalar density 
perturbations. We compute the angular power spectrum of curl-type polarization arising from small 
Faraday rotation due to a weak stochastic primordial magnetic field with a power-law power spec- 
■ trum. The induced polarization power spectrum peaks at arcminute angular scales. Faraday rotation 

is one of the few cosmological sources of curl-type polarization, along with primordial tensor pertur- 
bations, gravitational lensing, and the vector and tensor perturbations induced by magnetic fields; 
the Faraday rotation signal peaks on significantly smaller angular scales than any of these, with a 
power spectrum amplitude which can be comparable to that from gravitational lensing. Prospects 
for detection are briefly discussed. 
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A primordial cosmological seed magnetic field has been proposed to explain the existence of observed large-scale 
(~ 10 kpc), /iG strength, ordered magnetic fields observed in galaxies and galaxy clusters Q, 0- Cosmological 
perturbations induced by such a magnetic field and the corresponding cosmic microwave background temperature 
43 ' and polarization anisotropies have been the subject of a number of recent studies 0, 0, 0, 0, 0, H, El El El El 
Q-f El El El E] • The strength of such a hypothetical cosmological magnetic field can be significantly constrained by 
q • comparing model predictions and microwave background observational data • 

Several studies of the evolution of cosmological magnetic fields with power-law power spectra have been done E1EJ- 
The strongest constraints result from amplification of primordial magnetic fields via an inverse cascade mechanism, 
which takes power on small scales and transfers it to larger scales; this process is effective for fields with significantly 
\ increasing power on small scales 0, ■ Such "blue" power spectra result from magnetic fields generated in small 
scales at late time, for example, during cosmological phase transitions. In this work we consider general power-law 
magnetic fields, including fields generated during an early inflation-epoch, which are not strongly constrained by their 
subsequent evolution. Observational and theoretical constraints on primordial fields are discussed in more detail in 
" ' the concluding section of this paper. 

A distinctive method of constraining a cosmological magnetic field is to study the rotation of the microwave 
background polarization orientation due to the Faraday effect, as the radiation propagates to us in the cosmological 
magnetic field. A homogeneous magnetic field of strength 10~ 9 gauss induces a measurable rotation of order 1° at a 
frequency of 30 GHz j23j. Faraday rotation is the classic method of probing large-scale magnetic fields in the universe, 
using quasars as sources of polarized radiation, e.g. |24j . Faraday rotation of radiation from discrete point sources can 
constrain the primordial magnetic field power spectrum |25j . and can also probe the evolution of the magnetic field 
with redshift [2||. The microwave background in principle offers a great advantage because it provides a source of 
polarized light from a fixed redshift which covers the entire sky; with sufficiently precise measurements the projected 
magnetic field can be mapped in detail. The drawback is that the polarization of the microwave background is very 
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small, only a part in a million, and the Faraday rotation expected from a cosmological magnetic field is also small 
and so challenging to detect. 

In general, a polarization field has two independent components, known in the microwave background context as 
G and G [27j (for "gradient" and "curl"), or equivalently as E and B modes j2^|. It is well known that primordial 
scalar density perturbations induce only G-polarization, and that G polarization must arise from other production 
processes, notably primordial tensor perturbations |29l Isfij . gravitational lensing of the G-polarization component 
|3l| . and vector and tensor perturbations from, e.g., magnetic fields |l3j| . Faraday rotation, which rotates polarization 
orientations, can also be a source of G-polarization given an intial G-polarization field. 

In this paper we compute the power spectrum of microwave background G polarization resulting from Faraday 
rotation by a stochastic primordial magnetic field with a given power spectrum and helicity spectrum. We find 
that the power spectrum distinctively peaks on arcminute angular scales, significantly smaller than other commonly 
considered sources of G polarization. Previous work on the polarization power spectrum has considered the simpler 
but less realistic special case of a constant magnetic field (23, H^, |3l| , or considered a stochastic magnetic field but 
only calculated the power spectrum of the rotation measure rather than of the microwave background polarization 
|34| . Sophisticated analysis techniques for extracting microwave background power spectra from maps of the sky have 
been developed, so measuring the power spectrum of Faraday rotation and then confirming its scaling with frequency 
will be much more efficient than requiring frequency information from the outset to measure the rotation. Measuring 
the power spectrum of Faraday rotation also sidesteps any systematic errors associated with comparing observations 
at different frequencies; this is a significant advantage since these measurements will all be dominated by systematic 
errors. 

The rest of the paper is organized as follows. The next Section defines the primordial magnetic field power spectrum 
and helicity spectrum and explains the approximation we employ to separate polarization generation from polarization 
orientation rotation at the surface of last scattering. Section III derives the power spectrum of the Faraday rotation 
measure, which distills the Faraday effect of the tangled primordial magnetic field into a simple scalar function on 
the sky. We give an explicit demonstration that any helical component of the magnetic field does not contribute to 
Faraday rotation. Then given this rotation field and a primordial G-polarization field, Sec. IV derives an analytic 
expression for resulting power spectrum of G polarization. Section V numerically evaluates these expressions and 
presents the power spectra for both the rotation field and the G polarization for a range of magnetic field power 
spectra, and the final Section discusses these results in the context of future experiments to measure small-scale 
microwave background polarization. Two short appendices contains some useful but technical mathematical results. 



II. MAGNETIC FIELD MODEL 



A. Magnetic Field Power Spectrum 



We will assume the existence of a cosmological magnetic field generated during or prior to the early radiation- 
dominated epoch, with the energy density of the magnetic field a first-order perturbation to the standard Friedmann- 
Lemaitre-Robertson- Walker homogeneous cosmological spacetime model. The conductivity of the primordial plasma 
is high so we can work in the infinite-conductivity limit for all scales larger than the damping scale set by photon and 
neutrino diffusion. This results in a "frozen-in" magnetic field and a corresponding electric field E = — v x B where 
v is the plasma fluid velocity. Since the fluid velocity is always small in the early universe on the scales considered, 
the electric field is always negligible compared to the magnetic field and will not be considered further. Neglecting 
fluid back-reaction onto the magnetic field, the spatial and temporal dependence of the magnetic field separates, with 
magnetic flux conservation fixing the temporal dependence to be the simple scaling B(x, r)) — B(x)/a 2 where a is the 
scale factor and r\ conformal time. 

We also assume the magnetic field is a Gaussian random field. Taking into account the possible helicity of the field 
|35| , the magnetic field spectrum in wavenumber space is |l4l Il6| 

(B*(k)B,(k')> = (2^) 3 <5( 3 )(k- k'm.PBik) +ie ijl k l P H (k)], (1) 

where Pe(fc) and Pa{k) are the symmetric and helical parts of the magnetic field power spectrum, e^/ is the antisym- 
metric tensor, and the plane projector = 5{j — kikj satisfies PijPjk = Pik and Pijkj = with unit wavenumber 
components fcj = ki/k. We use the convention 




when Fourier transforming between real and wavenumber spaces. 
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The power spectrum P B (k) is related to the energy density of the magnetic field, while the helicity part P#(fc) is 
related to the spatial average (B • (V x B)) 14]. Transforming from an orthonormal basis {ei,e2,e3 = k} to the 
helicity basis [36| 

e ± (k) = --^=(ei ± ie 2 ), e 3 = k, (3) 
the power spectra P B (k) and P B {k) can be expressed in terms of magnetic field components in the helicity basis as 



the po-s 

Mm 



(B + (k)B+(-k') + B-(k)B-(-k')) = -(2Tr) 3 P B (k)6 (3 Xk-k') , (4) 
(fl+(k)fl + (-k / )-fl-(k)fl-(-k')) = (27r) 3 P H (fc)5 (3) (k-k'). (5) 

We describe both the symmetric and helical parts by simple power laws 

P B (k) = A B k nB , P H (k) = A H k n ». (6) 

If the magnetic field is generated in small scales, after inflation, the power spectrum index is constrained to be n_e > 2 
[37|. Also the amplitudes are generically constrained by P B (k) > \P B (k)\ ^HUES]], which implies n B > n B [l6L 
a field cannot support a fixed helicity in the limit of zero field strength. Note that some authors define spectral indices 
which correspond to our n B + 3 and Uh + 3 0, HU . 

As a phenomenological normalization of the magnetic field, we smooth the field on a comoving scale A by convolving 
with a Gaussian smoothing kernel / A (x) = N exp(— x 2 /2X 2 ) to obtain the smoothed field B A (x), and introduce average 
values of energy density 13] 

B x 2 ee (B A (x) • B A (x)) = — / dkk 2 e- k2x ' 2 P B (k), (7) 

n Jo 

and helicity (Itsf 

\ />oo 

i7 A 2 = A|(B A (x).{VxB A (x)})| = — / dkk 3 e- k2x2 \P H (k)\. (8) 

77 Jo 

Equations Q and © may be used to re-express the power-law power spectra of Eqs. © in closed form as [lli . IliL Il6| 

P H ( k ) = ±fIl _^ k<k D , (10) 

w 2 r (n H /2 + 2) kl H+3 

where the smoothing wavenumber k\ = 2tt/X. We make the approximation that both spectra vanish for all wavenum- 
bers k larger than a damping wavenumber kjy. We assume that the magnetic field damping is due to Alfven wave 
damping from photon viscosity, and the cut-off wavenumber to be 0, 0, Il3| 



n B +5 / d \ -2 / ; \ Ub+3 

2-9 xW^l J^r) h, (11) 



Mpc^v ' vio- 9 g; Vmpc" 1 

which will always be a much smaller scale than the Silk damping scale (thickness of the last scattering surface) for 
standard cosmological models (h is the Hubble constant in units of 100 km s _1 Mpc -1 ). 



B. Faraday Rotation 

A cosmological magnetic field at the epoch of last scattering will rotate the CMB polarization orientation in a 
given sky direction due to Faraday rotation (see, e.g .. [3 9llh For a review of CMB polarization theory see for 
computational methods and statistics, see |27l 123 . IllL |42 | . The time derivative of the orientation angle a of linearly 
polarized monochromatic radiation passing through a plasma in the presence of a magnetic field B(x, 77) is |23| 

da q 3 x e n e a 

c B (x, n, rj) = a ee - = ^—^ B(x, rj) ■ n. (12) 
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Here n is the propagation direction of the radiation and v the radiation frequency. The electron charge and mass are 
q and m e , and n e and x e are the total electron number density and ionization fraction. We normalize the scale factor 
a by setting ao = 1 today. We employ cosmological units with h = 1 = c. 

In terms of the comoving magnetic field B(x) = B(x, r])a 2 , Eq. (|12[) can be rewritten as 

3 

cjb(x, n, 77) = 2 f (x) B(x) ■ n (13) 

[w) u Q q 

with comoving frequency of the observed radiation Vq and the differential optical depth f = x e n e aTa. Here ctt — 
8TTa EM /(3m 2 ) is the Thompson scattering cross section with the fine-structure constant olem — q 2 ~ 1/137. We 
neglect inhomogeneities in the free electron density, assuming that f (77) is independent of x; then the rotation angle 
observed today is 

3 r m 

a(n, 770) = ia \2 2 / df/T(?/)B(x)-n. (14) 

Throughout the rest of this paper, B represents the comoving value of the magnetic field. 

Faraday rotation of the microwave background is a somewhat subtle problem, because the polarization is generated 
and rotated simultaneously in the region of the last scattering surface. In a rigorous treatment, these two effects 
must be computed together via the complete radiative transfer equations describing the evolution of polarization 
fluctuations. However, as long as the total rotation is small compared to tt/2, the total rotated polarization angle 
can be expressed simply as an average of the rotated polarization angle from each infinitesimal piece of path length 
through the surface of last scattering, neglecting depolarization effects. We make the simplifying approximation 
that any magnetic field component with a wavelength shorter than the thickness of the surface of last scattering is 
neglected. For these components, the rotation of polarization generated at different optical depths will tend to cancel, 
leaving little net rotation. (This assumption breaks down if the magnetic field strength is dominated by components 
on small scales. In this case, a more realistic power-law dropoff in Pb(/c) should be used instead of a sharp cutoff.) 
This is equivalent to imposing a cutoff scale on the magnetic field at the photon damping scale at last scattering. 
Then we can treat the magnetic field as constant throughout the rotation region, so that the total rotation, which is 
the sum over the rotations of each infinitesimal piece of generated polarization, can be expressed as the total rotation 
incurred by the polarization generated at some particular effective optical depth, which we denote by ?y r ot- 

Assuming that r) Iot is the conformal time corresponding to r = 1 (the actual value is a little lower since the 
polarization visibility function peaks at a lower redshift than the temperature visibility function), we can derive a 
simple expression for the approximate Faraday rotation in a given observation direction. Writing x = n(?7o — v) (i-e., 
putting the observer at the origin of the coordinate system) and Fourier expanding the magnetic field, Eq. (|14(l can 
be written as 

a(n ' %) ~ M2^vfq J d3fcB(k) ■ ne " k ' nA "' ( 15 ) 

where Ar/ = 770 — ry ro t- These approximations have split the polarization generation and rotation problems; we can 
treat the full problem as equivalent to the generation of polarization in the usual way, followed by rotation using 
an effective rotation screen just prior to the radiation reaching the observer. The rest of this paper is devoted to 
computing the properties of the effective rotation screen and the resulting power spectrum of microwave background 
polarization. 

III. ROTATION POWER SPECTRUM 

It is conventional to introduce a wavelength-independent measure of the rotation of the polarization orientation, the 
rotation measure R(n) = a(n)^Q. The two-point correlation function of the effective rotation measure for a stochastic 
magnetic field based on the approximate solution in Eq. (|15l) is 

(R(n)R(n')) ~ J d 3 k J d'k'^^e^' A "n^ (B*(k)B 3 (k')) (16) 

where A 77 = 770 - ?7rot ~ Vo- 

Prior to the ensemble averaging, it is convenient to decompose the vector plane wave into vector spherical harmonics 

M 

00 / 1 

B(k)e^ = £ E E A lm ( k ) Y im ( n ) (17) 

1=0 m=— i A=-l 
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where A labels the three orthonormal vector spherical harmonics defined by 

<><») = 7=1- v„ wm(11 ), 

YLV) = ny im (n). (18) 

Here V n is the two-dimensional covariant derivative orthogonal to the unit vector n, so the harmonics A = 1 and 
A = are transverse to n, while the A = — 1 harmonic is parallel to n. When the expansion l)17|l is substituted into 
Eq. 116(1 . only the A = — 1 terms contribute since the magnetic fields only appear through the contraction n • B. 
Explicit forms for the coefficients (k) are given in Ref. [3(| ; the one needed here is 



.4 



(-i) 
im 



in 



21 + 



jJ- 1 (vf+T) bwN) + ii-iN)]B(k) • 

- [(/ + l)j i+1 (fcr 7o )-/j ; _ 1 (fc7 ?0 )]B(k).Y i ( m 1) *(k)}. (19) 
The second term proportional to Yt 1 ' is zero since the magnetic field is divergenceless, k • B(k) = 0, leaving 

4? = ^ V^TT)^B(k) • Y«*(k). (20) 
With these definitions the rotation angle becomes 

Q(n ' ??o) " 32^ / AE'^^^Bfk) • Y«*(k)Y im (n), (21) 



and for the magnetic field spectrum of Eq. Q we obtain 
(R(n)R(n')) ~ J dkk 2 P B (k) J dfi fc 



x E E *''~V*(* + W + l) J ' (fc g J ;l fc??o) ^(n)y, m ,(nO Y|:»(k) ■ Y&T(fc). (22) 

Zm Cm' 

In this result, the term proportional to fcjfej in the projector Pjj in the magnetic field power spectrum in Eq. 

vanishes because k • YS^ = 0. The term proportional to the magnetic field helicity spectrum Pji(k) is also identically 

zero because [k x Y,^ (k)j • Y^j^k) — Y^(k) • Y^^k) = 0. This result is in agreement with Ensslin and Vogt 01 
and Campanelli et al. |3J| , who conclude that Faraday rotation cannot be used to reconstruct the helical part of the 
magnetic field spectrum. The opposite conclusion of Ref. |l4j is erroneous. Physically, in real space, non-zero helicity 
of the magnetic field only affects the off-diagonal elements in the correlation matrix of the magnetic field components, 
while the Faraday rotation is due only to the diagonal component corresponding to the propagation direction; see |43j 
for a more detailed discussion. 

Using the orthogonality of vector spherical harmonics, 

J c%Y W * (k) • Y^, (k) = Sxx, 5w <W , (23) 

and the usual spherical harmonic summation formula 

An 



P,(n • n') = — ]T Y^(n)H m (n'), (24) 

771 — — I 

Eq. J23) simplifies to 

(i?(n)i?(n,)> * i2W^ + 1)p,{n ■ n ° / dkk2pB{k) (^) 2 ' (25) 
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with corresponding multipole moments 

91(1 + 1) B\ /A^" B+3 r™> 

~ (4tt)V r(n s /2 + 3/2) {r,o J ./„ 

where xjj — kjjrjo and the multipole moments are defined via 



^Swd^lrl / (26) 



<i2(n)fl(n')> = £ -^-Cfi^n • n'). (27) 

The rotation angle power spectrum is simply given by the rescaling 

Cf = v^Cf. (28) 



Note that Refs. [Fa. l4q have incorrect prefactors in their corresponding expressions for the rotation multipoles C/*, 
due to the non self-consistent choice of the units. 

This rotation multipole expression contains a sharp short- wavelength cutoff kjj ; in reality, the effective cutoff will be 
smoothly spread over a range of scales. To prevent unphysical oscillations in the integral of Eq. (|26(l , and to simplify 
the numerical evaluation of the integral, we replace the oscillatory function jf(x) by half of its envelope, l/(2x 2 ), for 
all x larger than the second zero of ji(x). Then the tail of the integrand is just a trivial power law, and the total 
integral will not show oscillations, as would be expected from a more realistic cutoff function. 

IV. C-POLARIZATION FROM FARADAY ROTATION 

The preceding Section gives an expression for the power spectrum of an approximation to the polarization orientation 
Faraday rotation field a(n) which can be applied to the unrotated polarization field of the microwave background 
radiation. If the unrotated polarization field arises only from scalar perturbations it will only have a nonzero G 
polarization component. Faraday rotation will induce non-zero C polarization; here we compute the C-polarization 
power spectrum. We use the differential geometry formalism of Ref. |27| for this computation, and work with tensors 
defined on the two-dimensional spherical manifold representing the sky. 

We represent the unrotated polarized CMB with the tensor field P ab (h) on the two-dimensional sphere orthogonal 
to the direction vector n, and decompose this in terms of tensor spherical harmonics in the usual way 



P - \n G V G Jun C V C 

ab ~ \ U lm (lm)ab "T" "'Im 1 (lm)ab 



(29) 



Here the usual Im indices are enclosed in parentheses to distinguish them from the ab tensor indices. Explicit forms 
for the orthonormal tensor spherical harmonics are 

Y ( f m)ab (n) = N t (Y (lm y Mb (n) + ^g ab l(l + l)Y {lm) {n)^ , 

Y$ m)ab {n) = y (Y {lm y. ac (n)e c b + Y {lm):bc (n)e c a ) , (30) 

with the normalization factor Ni = (2(1 — 2)!/(Z + 2)!) 1 / 2 . Here g a b is the metric on the sphere, e a b is the antisymmetric 
Levi-Civita tensor, and the covariant derivative on the sphere is A(n) . n = V n ^4(n), using the notation of the previous 
Section. For polarization from scalar perturbations, af m = |29t l30j|. Primordial tensor perturbations, or the vector 
and tensor perturbations induced by magnetic fields, will give nonzero af m . 

A rotation of the polarization orientation by an angle a is represented by the transformation 

P' ab = R a c P cb (31) 

with the rotation operator given by 

R a c - cos(2a). 9a c + sin(2a)e a c . (32) 

In the limit of small rotation angle, a«l, the rotation operator reduces to R a c ~ g a c + 2ae a c , which is linear in a. 
In our case the rotation field varies with direction a = a(n), so we decompose it into spherical harmonics 

«( n ) =^2b lm Y (lm) (n). (33) 
Im 



7 



We now use Eq. I|29(l to expand the resulting polarization field in tensor spherical harmonics. Inverting for the 
coefficients of the rotated C polarization we have, 



7 C ' 
''Irn 



and likewise for the G-polarization multipole coeffecients 



J dn n P' ab (n)Y^*(n) = 2 J dQ n a(n)e a c P cb (n)Y^f*(n) 

2 E b h mi a? 2m2 [ dQ a Y (limi) (n) e a c Y ( % m2) cb (n)Y^f *(n), (34) 



a frn = 2 E Yl bl ^ a Zm 2 J ^n% mi )( n ) € a 



C Y ( Z m2) Jn)Y^*(n). (35) 



The integrals over three spherical harmonics can be evaluated through repeated integrations by parts, after making 
use of the form of the tensor spherical harmonics and the identities 

tactbd = gabged — g a d,gcb, Cca^b = Bab = —eacC C bi e ab-c = (36) 

which eliminates the product of antisymmetric tensors in the integrand. We also need the eigenvalue equation 
Y(im):a' a ( n ) = — '0 + i )^ / (Zm)( n ) an d use the integration by parts formula 

dn n A :a B = - J dn n AB.. a , (37) 

which has no surface term contribution as long as the integration is over the entire sphere. Using the explicit forms 
for Y^ rn)ab {n) and Y^ m)ab (n) in Eqs. l|3"0"j) and using Eq. (|3"6*|l results in 

J dn n Y (hmi) (n) e a c Y ( f 2m2)cb (n)Y^ } b *(n) = ^ J dfi„Y (Iim0 (n)Y ( , 2ma):6 :6 (n)Y ( ; m):a :0 (ii) 

-N l2 m J ^ n r (iimi) (n)Y( i2m2):ah (n)r ( t m) :Qb (n). (38) 

For G <-> C in the integrand the result is the same except for an overall sign change. The derivatives in the first 
integral just add 1(1 + 1) factors. 

The second integral is more involved, and is computed in Appendix A. Substituting Eqs. (| AT|) and IjAlJ) into Eq. I|38(l 
and then into Eq. (|34J) gives finally 

a?m = ^EE N l2 K(l,h,l 2 )b hmi af 2m2 J dnY (hmi) Y (l2m2) Y( lm) , (39) 

limi l 2 rri2 

where 

K (I, Ji, f 2 ) = ~\ {L 2 + Lj + L 2 2 - 2L X L 2 - 2^L + 2Li - 2L 2 - 2L) , (40) 

with L = 1(1 + 1), L 1 = l\(li + 1), and L 2 — l 2 (l 2 + 1). The integral over three spherical harmonics is a well-known 
expression in terms of Clebsch-Gordan coefficients, Eq. IjAlfl . The coefficient a^' is given by the same formula with 
the replacements af 2Tn2 -» af 2 77l2 and Ni -> -iVj. 
To derive the power spectrum, we need to evaluate 



,c l*„C i\ 



Cftu'Smm', (41) 



where the definition of Cf holds for any statistically isotropic distribution; here we assume that all relevant cosmo- 
logical quantities satisfy this condition. Assuming that the rotation field and the original temperature polarization 
field are uncorrelated, we have 

( b l 1 m 1 b hm2 a l 3 m3 a l 4 m i ) = &h h 5 mi m 2 ^hU ^rn 3 m 4 C;" C; 3 • (42) 

where Cf is the angular power spectrum of the rotation field and CF is the G-polarization power spectrum of the 
unrotated polarization field just after decoupling. Using the expressions given in Eqs. (|39|l and (|42|) in Eq. I|41|) results 
in a long expression with sums over (pi and l 2 m 2 . This can be further simplified using the Clebsch-Gordan identity 

E Ci™ ni i 2m2 C l ll ™ ll2 „ l2 = dw 5mm', (43) 
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FIG. 1: The angular power spectrum of the Faraday rotation angle, Cf, induced by a stochastic magnetic field. The curves in 
order of decreasing amplitude on the right side of the plot correspond to magnetic field power spectral indices tib = 2, 1, 0, 
— 1, and —2. The magnetic fields have been normalized to a nanogauss at the smoothing scale A = 1 Mpc. 



and we find 

C? = Nf £ N?K(l, luhfCgCZ ^(2/+ i) + ^ (^o) 3 ■ ( 44 ) 

This is our basic result for the C-polarization power spectrum induced by Faraday rotation of the primordial G- 
polarization; we numerically evaluate it in the following Section. (Numerical techniques for evaluating the remaining 
Clebsch-Gordan coefficients are given in Appendix B.) Note that the more familiar multipole coefficients Cf and Cf 
psl ] are just a factor of two larger than Cp and Cf , respectively. 

We note that, for a stochastic magnetic field, the cross-correlations between the intrinsic temperature and G- 
polarization fluctuations from non-magnetic perturbations and the C-polarization from Faraday rotation result in 
Cff = 0. These multipoles are linear in the magnetic field so the ensemble average vanishes. Non-zero correlations 
result only from the intrinsic T fluctuations arising from the magnetic field itself and the C-polarization from Faraday 
rotation, but this case is not considered here. The Cp c correlation will also be zero unless the magnetic field has 
non-zero helicity, which follows from the invariance of (QU) under Faraday rotation by a stochastic magnetic field 

El 1 . 



non-zero off diagonal correlation will apppear in the case of homogeneous magnetic field Ref. l33l 
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FIG. 2: The C-polarization power spectrum of the microwave background induced by the Faraday rotation field in Fig.0 again 
with the magnetic field normalization scale A = 1 Mpc. 



NUMERICAL RESULTS 



The rotation angle power spectrum, Eq. I)28[l. is shown in Fig. ^for values of the magnetic field spectral index 
ranging from tlb = —2 to fig = 2. The plotted rotation power spectra are for a magnetic field with B\ = 10~ 9 G 
and an observation frequency of v = 100 GHz; the rotation power scales like B\v~ A . We have assumed a cutoff scale 
fcrj = 2.0 Mpc -1 approximately corresponding to the Silk damping scale. The polarization has a root-mean-square 
rotation angle given by 



/ 2\V2 



E 



21 + 1 

4-7T 



C? 



1/2 



(45) 



At v — 100 GHz and B\ = 10 -9 G with A = 1 Mpc, a w 0.3° for all considered values of ns- This is consistent with 
a found in Ref. |23j for a constant magnetic field of strength B\ (around 1.6° at v = 30 GHz). For negative values 
of ns + 1, the magnetic field power spectrum grows with length scale, so the largest rotations are seen at large scales 
corresponding to small I values. The opposite is true for positive values of Ub + 1, which increase in power at smaller 
scales and larger / values. 

The power spectrum of C-polarization induced by this Faraday rotation power spectrum is displayed in Fig. [21 
This power spectrum is also for B\ — 10~ 9 G and v — 100 GHz, scaling like B\v~ A . The most notable feature of 
the C-polarization spectrum is its peak at small angular scales at I > 10 4 . The cutoff at I = 16000 is sharper than 
in a more realistic damping model where it occurs over a range of scales somewhat past fc^, but the peak position 
and amplitude would be largely unchanged by a different damping prescription. We also assume a standard ACDM 
cosmological model and compute the Cf power spectrum appearing in Eq. il4*4")l to I = 5000 using CMBFAST p5|. 

The peak angular scale is readily understood: at the last scattering surface Faraday rotation induces polarization 
fluctuations on angular scales corresponding to the characteristic wavelengths of the stochastic magnetic field; just sub- 
sequent to recombination the polarization is imprinted with fluctuations on scales much smaller than the quadrupole 
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scale which induced the polarization fluctuations in the first place. Free streaming then shifts these fluctuations to 
even smaller scales. So the Faraday rotation power spectrum provides a distinctive signature of primordial magnetic 
fields. 

This computation only applies to small rotation angles, which for a given magnetic field translates into a limiting 
lower observation frequency. The root-mean-square rotation a, Eq. Q45[l. scales like B\Vq 2 . For a field strength of 
10~ 9 G, observations at v = 30 GHz yield a mean rotation of around 3°, well within the small rotation limit, but 
observations at v = 10 GHz would give mean rotation angles of around 30°, for which the small- angle approximation 
employed in Eq. (|34|l clearly is not valid. If the rotation angle becomes large, depolarization effects become important 
and the mean polarization amplitude is reduced |46| . 

The other sources of cosmological C-polarization have been widely discussed. The first is tensor (or vector) primor- 
dial perturbations. In particular, an epoch of inflation in the early universe necessarily produces nearly scale-invariant 
tensor as well as scalar perturbations, and for many inflation models these tensor modes are observable through the C- 
polarization they produce 01 • This class of tensor perturbation signals all have the property that their C-polarization 
power peaks at large scales, / < 100. The other generic source of C polarization is gravitational lensing of the pri- 
mordial G polarization from scalar perturbations j^jj, which peaks at angular scales in the region of I « 1000. The C 
polarization from Faraday rotation peaks at significantly smaller angular scales, well separated from either the tensor 
or lensing signals, and would thus be easy to separate. 



VI. DISCUSSION AND CONCLUSIONS 



Primordial magnetic fields produce C-polarization fluctuations directly via the vector and tensor perturbations they 
induce in the primordial plasma |l3j . For a field with B\ = 10 -9 G, the peak polarization amplitude is, for example, 
l 2 Cf « 10~ 13 at I = 1000 for n = -2 [l]J. The Faraday rotation signal from the same magnetic field gives a peak 
polarization amplitude of around 10" 14 for v = 30 GHz and (neglecting the depolarization effect at significant rotation 
angles) 10~ 12 for v = 10 GHz. The direct polarization fluctuations peak at significantly larger scales than those from 
Faraday rotation: primordial magnetic fields have a double signature, with two polarization peaks. In comparison, 
the G polarization power spectrum from primordial density perturbations peaks at a comparable amplitude of around 
Z 2 Cp w 3 x 10~ 12 at an angular scale of around I — 1000, while the C polarization from gravitational lensing has an 
amplitude on the order of l 2 Cf « 10~ 14 at I — 1000. The C-polarization from Faraday rotation for realistic primordial 
magnetic fields is small but potentially measurable, in the same ballpark as detecting the lensing polarization signal. 

Several claims have been made that primordial magnetic fields of the amplitudes considered here are ruled out. 
These claims are not all compelling; we review some here. Caprini and Durrer have asserted remarkably stringent 
limits on primordial magnetic fields through their conversion into gravitational waves, which then are constrained 
by the expansion rate of the universe at nucleosynthesis 48] . But the expansion rate of the universe is the same 
whether energy density is converted from magnetic fields to gravitational waves or not, since the energy density of 
both scale the same way with the expansion of the universe. So the actual constraint is on the total radiation energy 
density in the magnetic field, which is constrained to be about 1% of the total energy density in the usual manner 
that nucleosynthesis limits extra neutrino species. The corresponding limit on the total comoving mean magnetic 
field strength is around 10 -8 gauss, not the 10 -27 gauss claimed in [48fl . 

Potentially more interesting limits come from magnetic fields in galaxy clusters. It is believed that clusters contain 
magnetic fields with u,G field strengths, from Faraday rotation of background polarized radio sources observed through 
the cluster (e.g., |43. l50|'). Vogt and Ensslin [5l| have recently analyzed several cluster rotation maps, concluding 
that fields of at least a few fiG are present, with a lower limit on the steepness of the magnetic field power spectrum 
in the clusters. On the other hand, simulations starting from constant comoving magnetic fields of 10~ n G show 
clusters generating fields sufficiently large to explain Faraday rotation measurements |20ll52j |. Banerjee and Jedamzik 
|lj| have claimed that these results rule out primordial fields with comoving amplitudes as large as 10 -9 G, arguing 
that they would overproduce fields in galaxy clusters today. We feel this is a potentially powerful argument, but 
premature. First, simulations have shown that primordial fields which are 1% of those considered in this paper may 
lead to observed cluster fields, but magnetic field generation is a non-linear process and it is not obvious that the final 
cluster fields will scale linearly with the initial fields; we know they will definitely not as equipartition field strength 
is approached in the cluster. While it seems likely that cluster fields are generally below equipartition value 53j, 
simulations over a wide range of initial field strengths and configurations must be performed to rule out certain initial 
field configurations. Second, we do not fully understand galaxy clusters, and it is not immediately clear whether the 
cluster simulations to date include all of the relevant physics. Third, the Faraday rotation signal in galaxy clusters 
depends on not only the field strength but also its configuration: it is conceivable that larger initial fields might lead 
to larger cluster fields, but that much of the field strength would be concentrated at small scales due to more efficient 
turbulent cascades, and small-scale fields contribute much less to observed Faraday rotation measurements. It is likely 
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that further cluster simulations will resolve these issues; until then it is important to examine other techniques for 
probing magnetic fields like the one outlined in this paper. 

Another constraint on primordial fields comes from field dissipation on small scales inducing spectral distortions 
in the microwave background radiation. Field strengths greater than 3 x 10~ 8 G on scales of 400 pc are ruled out. 
This limit constrains primordial magnetic fields which are blue, n > 0, with most field strength on small scales. For 
the fields considered in this paper, on scales larger than 1 Mpc, such constraints are important if the magnetic field's 
power spectrum remains a power law over a range of 10 4 in length scale. As this may not be the case, we present 
results in this paper for various values of n, but the reader should note that pure power-law fields with n > 0.5 must 
have amplitudes smaller than 10~ 9 G on Mpc scales. 

A signal we have not computed here explicitly, but which follows simply from the formulas derived above, is 
the G-polarization signal arising from the Faraday rotation of a primordial C-polarization signal. Of course, scalar 
perturbations provide no such primordial signal, but any magnetic field does [l3| . It is likely that the magnetic field 
polarization signal is small compared to the G-polarization from primordial density perturbations, but the rotated 
G polarization appears, like the G polarization, at small angular scales of / 10000. At these small scales, any 
direct primordial signals will be suppressed due to diffusion damping, so this G polarization signal would provide 
confirmation of the corresponding C polarization rotation signal in Fig. [21 

Additionally, comparison of microwave background imprints from Faraday rotation and intrinsic fluctuations can 
probe the helicity of the primordial field. Here we have demonstrated explicitly that any helical part of the magnetic 
field gives zero Faraday rotation. Helical stochastic magnetic fields do, however, produce non-zero temperature and 
polarization fluctuations directly 0, . By comparing the microwave background power spectra at scales of 20' 
with the polarization power spectra at scales of 1', the helical power spectrum Pjj (k) can be probed separately from 
the non- helical power spectrum Pg(fc). 

In this paper, we have modelled the Faraday rotation signal by separating the polarization generation and rotation 
processes through a simple approximation. A more accurate calculation demands a full numerical evolution of the 
coupled Boltzmann hierarchy of equations describing photons, electrons, and magnetic and gravitational fields. Such 
a code has been developed for the direct temperature and polarization fluctuations from magnetic fields and at 
first glance it would appear to be a simple matter to modify such a code to include Faraday rotation, including a more 
accurate treatment than presented here of combined polarization generation and rotation, and damping at small scales 
or large magnetic field strengths. However, currently used microwave background Boltzmann codes use polarization 
basis functions corresponding directly to G and C polarization modes. While this provides formal simplicity and 
conceptual clarity, Faraday rotation, described by the rotation operator Eq. i|32|) . does not correspond to a simple 
rotation of G polarization into C polarization, but is naturally expressed as a rotation mixing the Q and U Stokes 
parameters. A straightforward numerical implementation of Faraday rotation requires using a code employing the 
Stokes parameters as variables. Techniques for doing this and reconstructing the G and C power spectra are known 
(see, e.g., 01) an d were used in a number of older microwave background polarization codes. 

The polarization signals discussed in this paper have an amplitude of /iK or smaller, on angular scales of 1' or smaller. 
The combination of small amplitude and angular scale make detection of this potential signal challenging. For a given 
magnetic field strength, the Faraday rotation power spectrum scales like v with observing frequency, so going to a 
low enough frequency will compensate for smaller magnetic field strengths. But lower-frequency observations at the 
same angular resolution require larger experiments, and to obtain arcminute resolution at frequencies of 30 GHz or 
below over a reasonable field of view likely requires interferometric experiments (the corresponding diffraction limit 
would require a single-dish diameter of 50 meters). No polarized interferometers acting at these frequencies and 
ang ular scales are currently envisioned, but polarization experiments have only just begun to detect any signals at all 
|54l l55l l56l l57j , and we can expect rapid advances over the coming decade. The new generation of millimeter- wave 
frequency high-resolution experiments for temperature fluctuations (e.g., |5£j^ will help lay the technical basis for 
complementary polarization observations. 

The Faraday rotation signal in Fig- EI is & demanding experimental target, and its amplitude is highly uncertain 
(and could certainly be zero). But we will not be obliged to go directly from ignorance to probes of these small scales, 
because cosmological magnetic fields strong enough to induce detectable Faraday rotation will also cause temperature 
and polarization fluctuations at larger scales. If primordial magnetic fields are present, we will suspect their existence 
long before the Faraday rotation power spectrum is probed directly, and the existence of detectable Faraday rotation 
can be verified through multi-frequency polarization observations along individual lines of sight. If the coming five 
years strengthen current suspicions of a possible primordial magnetic field into more solid evidence, then the Faraday 
rotation signal calculated in this paper will become a compelling experimental opportunity to characterize a new and 
poorly understood relic of the early universe. 
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APPENDIX A: INTEGRALS INVOLVING THREE TENSOR SPHERICAL HARMONICS 

The evaluation of the needed integrals over three tensor spherical harmonics can be built on the integral over three 
usual (scalar) spherical harmonics using repeated integration by parts. The integral over three spherical harmonics is 
well known |3j| , 

AOV V V* - ( ( 2?1 + 1 )( 2 ^ 2 + plO n lm /A-n 

«" * (hm 1 )*{l 2 m 2 ) :r (lm) ~ I ^k{21 + 1) / U llOhO^ lirnihrn? > 

where C) rr L , „, is a Clebsch-Gordan coefficient. 

Now we need to compute the second integral in Eq. H38fl . 

h = J rfr» n r (/imi) (n)r (/2m2):ah (n)y ( t m) :Qb (n). (A2) 

First consider the simpler integral 

h = |^ n r (iimi) (n)y (;2m2):a (n)r ( t 3m3) :Q (n). (A3) 

This can be evaluated by integrating by parts three times: first on the derivative of Yn 3fn3 -\, which gives two terms, 
for one of which the derivatives change into h(h + 1)- Now for the other term, integrate by parts a second time on 
the derivative of y| l!m2 ), and then follow by a third integration by parts on the derivative of Yn imi \. The integral 
which is left is the negative of the original integral; solving for this integral results in 

h = ±(-L 1+ L 2 + L 3 ) J dn n Y (limi) (n)Y ihm2) (n)Y(i sm3) (n) (A4) 

with the abbreviation Li = + 1). The remaining integral is just Eq. 1A1(I . 

Using Eq. I|A4|I , the integral l|3"5|) can be evaluated using the same general strategy, with three integrations by parts. 
In this case, we also need to compute Y a b' b = Yb a b . To simplify this, we need to commute the two outer derivatives. 
The curvature tensor for this 2-dimensional maximally symmetric space is just 

Rabcd = 9ac9bd ~ 9ad9bc- (A5) 

Thus we have 

Y ab b = Y :ba - b = Y-- b .. ba + R dba b Y- d = -[1(1 + 1) - l]Y a . (A6) 



Now we have all the pieces needed to evaluate Eq. (|A2(1 : after three integrations by parts, we are again left with the 
original integral, plus surface terms, like in the case of Eq. I)A3|I . The surface terms can all be evaluated in closed 
form with the help of Eqs. l|A6j) and (|A4|I . giving 

h = \{L 2 +L-L t - 2)(L 2 + L-L{)j <m n Y {hmi) Y {hm2) Y { * m) . (A7) 

Again, the remaining integral is given by Eq. IjAlfl . 



APPENDIX B: CLEBSCH-GORDAN COEFFICIENT NUMERICS 



The Clebsch-Gordan coefficients in Eq. (|44|l can be evaluated in closed form as |3f 



°a060 — 



(-l)9- C V2^+Tff! 

(g-a)Kg-b)\(g-c)\ 



(2 5 -2a)!(2 5 -26)!(2.g-2c)! 
(^+1)1 



1/2 



(Bl) 
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which is valid for a + b + c = 2g with g a positive integer; a, 6, and c must also satisfy the triangle inequalities. If 
a + b + c is odd or if a, 6, and c cannot form the sides of a (possibly degenerate) triangle, the coefficient vanishes. 
This coefficient can be directly evaluated numerically, as long as the various factors are taken in an order to prevent 
overflow and/or underflow errors. However, direct evaluation is not fast enough for efficient evaluation of Eq. I|44f) . 
since performing the sums up to / values of several thousand demands millions of Clebsch-Gordan coefficients, most 
with very large indices. But the two-term Stirling asymptotic expansion 

nl ~ n"e~"V2^ (l + — ] , n oo, (B2) 



12n, 

is accurate to 0.1% even for n — 1. Applying this expansion to all of the factorial factors in Eq. (|B1|) gives the 
approximation 

when g~a, g~b, and 5 — c are all non-zero. This approximation is accurate to just over 1% for the case a = b = c = 2, 
and better than 1% for all other needed terms. For the degenerate case a + b — c, 

<&>!.._!_ f 1+ f 1 _ 1 _ n (9o6) -v* . (B4) 



2c + 1 2/tt \ 2gJ \8g 8a 

along with cyclic permutations for the other two cases; this approximation is accurate to better than a part in 10~ 4 
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